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ABSTRACT 

We calculate Gamma-Ray Burst (GRB) afterglow light-curves from a relativistic jet of initial opening angle 9q, 
as seen by observers at a wide range of viewing angles, dohs, from the jet axis. We describe three increasingly 
more realistic models and compare the resulting light-curves. An observer at 6'obs < 6*0 should see a light curve 
very similar to that for an on-axis observer. An observer at 9ohs > do should see a rising light curve at early 
times, the flux peaking when the jet Lorentz factor ~ 1/6'obs- After this time the flux is not very different from 
that seen by an on-axis observer A strong linear polarization (<; 40%) may occur near the peak in the light curve, 
and slowly decay with time. We show that if GRB jets have a universal energy, then orphan afterglows associated 
with off-axis jets should be seen up to a constant dohs, therefore the detection rate of orphan afterglows would be 
proportional to the true GRB rate. We also discuss the proposed connection between supernova 1998bw and GRB 
980425. 

Subject headings: gamma rays: bursts — ISM: jets and outflows — radiation mechanisms: nonthermal 



1. INTRODUCTION 

Gamma-Ray Bursts (GRBs) are explosions which release 
roughly 10'^^ erg in the form of kinetic energy of highly rela- 
tivistic material (Frail et al. 2001, Panaitescu & Kumar 2001)^ 
Many GRBs appear to be highly non-spherical explosions, as 
evidenced by a nearly-achromatic break in the light-curve (e.g. 
Harrison et al. 1999; Stanek et al. 1999). Highly relativistic 
jets are "visible" when our line of sight is within the jet aper- 
ture (6'obs < (^o), otherwise, because of relativistic beaming of 
photons away from our line-of-sight, the object is too dim. As 
the jet decelerates, the relativistic beaming becomes less severe 
and the emission from the jet becomes detectable to observers 
at larger viewing angles. 

In this Letter we study the afterglow light-curves for off-axis 
locations (^obs > 0), focusing on observers lying outside of the 
initial jet opening angle (6'obs > 6'o). Granot et al. (2001) have 
shown that the light curve seen by an observer located within 
the initial jet aperture (6'obs < 6'o) is very similar to that for an 
on-axis observer (0obs = 0). Dalai et al. (2002) and Rossi et 
al. (2002) have presented simple models to calculate the flux 
in this case. We reanalyze these models in §2.1 and consider 
more realistic models in §2.2 & §2.3. Moderski, Sikora and 
Bulik (2000) have calculated off-axis light-curves with a more 
complex model, similar to that presented in §2.2. 

In §3 we calculate the temporal evolution of the linear po- 
larization for various 6'obs- In §4 we analyze the prospects of 
using the detection rate of orphan afterglows to estimate the 
collimation of GRB jets. In §5 we analyze the suggestion of 
Woosley, Eastman, & Schmidt (1999) that a relativistic jet em- 
anating from the SN explosion and pointing away from us could 
explain the observations. 

2. MODELING THE OFF AXIS EMISSION 



order to improve our understanding of the underlying physics 
and in order to check how general the results are, we explore 
three different models with an increasing level of complexity. 

2. 1 . Model 1: A Point Source at the Jet Axis 

We begin with a simple model, where for 6'obs — the 
light curve follows the results of simple jet models (Rhoads 
1999; Sari, Piran & Halpern 1999, hereafter R-SPH99), and 
for 0obs > the light curves are calculated assuming the emis- 
sion is from a point source that moves along the jet axis. The 
on-axis light curve exhibits a jet break at (R-SPH99): 

ijet = 6.2(1 + z)(S52/7io)'/'(^o/0.1)'/' hr , (1) 

where E52 is the isotropic equivalent energy in units of lO'^^ 
erg, no is the ambient density in cm^'^ and z is the cosmolog- 
ical redshift of the source. Att< tjct, -Fl/ (6'obs = 0) is taken 
from Sari, Piran and Narayan (1998), while alt > tjet the tem- 
poral scalings of the break frequencies and peak flux change 
according to R-SPH99. The observed flux density from a point 
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where L'^, and i^' are the spectral luminosity and frequency in 
the local rest frame of the jet, and are the angular and 
luminosity distances to the source, 7 = (1 — is the 

Lorentz factor of the source and 6 is the angle between the 
direction of motion of the source and the direction to the ob- 
server in the observer frame (in our case 6 = Oohs)- Since 
t/t' « dt/dt' = v' /v = (1 + 2)7(1 - 13 cose), where t and v 
are the observed time and frequency, we obtain that 



In this section we calculate the afterglow light curves of jet- 
ted GRBs, as seen by observers at different viewing angles, 
6'obs, w.rt the symmetry axis of the jet. For simplicity, we con- 
sider only a jet propagating into a homogeneous medium. In 

^ Most of the information we have about GRB explosions is only for the so-called long bursts, lasting more than a few seconds 
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where tg and 1^0 are the observed time and frequency for an 
observer at 0obs = ^- One therefore obtains that 

F,{dobs,t)^a''F,/a{0,at) , (4) 

where, for simpHcity, we take 7 = 6*0^^ I^o /^jct] '^^^ at to < tjct 

and 7 = 6io"^[io/<jct]"^/^ at io > ^jct- 

The light curves obtained using equation Q are shown by 
the dashed Unes in Figure |l[ At first 76'obs ^ 1 and « 
(7^obs)~^ is the dominant term in equation H, giving a sharp 
rise in the light curve. Once 7 becomes <^ 6^^^^ the flux begins 
to decay, asymptotically approaching the on-axis light curve. 
The light curve for off-axis observers peaks when 7^1 /6'obs- 
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Fig. 1. — B-band luminosity for models 1 (dashed lines) and 2 (solid 
lines), for 6*0 = 5°, 6»obs = (0, 1, 2, 3, 5)6lo, E52 = 80, no = 1, 
p = 2.5, €g = 0.01, Ee = 0.1, where eg (ee) is the fraction of the 
internal energy in the magnetic field (electrons) and p is the power law in- 
dex of the electron energy distribution. Model 1 is scaled down by a factor 
of 2.5, to help compare between the two models. 

The main advantage of this model is that it is very simple, and 
nevertheless gives reasonable results for 6'obs ^ 20o- Its main 
drawback is that it is not physical for 6'obs ^ at t ^ tj^t- 
This is because for ^obs < Oq the observed radiation is ini- 
tially dominated by emission from the material within an angle 
of 1/7 < 6'obs around the line of sight, while in model 1 the 
emission is always only from along the jet axis, and therefore 
the dominant contribution to the emission is missing, until the 
time when 7 ^ 1 /6'obs- This problem is overcome by our next 
model. 

2.2. Model 2: A Homogeneous Jet 

This model is described in Kumar & Panaitescu (2000), and 
here we briefly point out its main features. The Lorentz factor 
and energy density per unit solid angle are considered to be in- 
dependent of angle 6 within the jet aperture. The decrease of the 
Lorentz factor of the jet with time is calculated from the mass 
and energy conservation equations, and the sideway expansion 
speed of the jet is taken to be the local sound speed. 

The radiation calculation includes the synchrotron and 
inverse Compton processes, and the synchrotron spectrum 
is taken to be piece-wise power-law with the usual self- 
absorption, cooling and the synchrotron peak frequencies cal- 
culated from the electron spectrum, magnetic field strength and 
the radiative loss of energy for electrons. The observed flux is 
obtained by integrating the emissivity over equal arrival time 
surface (e.g. Kumar & Panaitescu 2000). 



The light curves of model 2 are shown by the solid lines in 
figure |l} The flux density in the decaying stage is slightly higher 
for larger viewing angles 6'obs- This effect occurs since at this 
late stage the whole jet is visible, and for larger 6'obs the radia- 
tion from a given radius arrives at the observer at a latter time, 
on average. Therefore, for a given observed time, larger 6'obs are 
dominated by emission from smaller radii, resulting in a larger 
flux density. At a few hundred days, the light curves begin to 
flatten due to the transition to the non-relativistic regime. 

The light curves for 6'obs <, Oq are very different from 
model 1 (and more realistic). Furthermore, the light curves for 
6'obs < 00 are very similar to 6'obs = in this model. Since the 
jet is homogeneous, the ratio of the observed flux for 0obs < 6*0 
and 6'obs = 0, may be approximated by the ratio of the areas 
within the jet, that are within an angle of I/7 around the di- 
rections to these two observers (which never decreases below 
1/2). 

We notice that the lights curve of model 1 for 6'obs /^o = 1,2 
are much closer to the light curves of model 2 for ^obs/^'o = 
2, 3, respectively, than to the light curves for the same view- 
ing angles. This is so because the emission for an observer 
outside the jet opening angle is dominated by the point in 
the jet closest to the observer Therefore, model 1 will be- 
come much more realistic (and just as simple) if one would use 
= max(0, 9ohs — (^a), rather than 9 = 9ohs, in equations ^and 

i 

The main advantage of model 2 is that it provides realistic 
light curves in a very reasonable computational time, making 
it very convenient for performing detailed fits to observations 
(e.g. Panaitescu & Kumar 2001). Its main drawback is a rel- 
atively simple treatment of the dynamics, which causes some 
differences in the light curves, compared to our next model. 

2.3. Model 3: 2D Hydrodynamical Simulation 

This model is described in Granot et al. (2001). The jet 
dynamics are determined by a 2D hydrodynamical simulation, 
with initial conditions of a wedge taken from the spherical self 
similar Blandford-McKee (1976) solution. The light curves for 
observers at different 9ohs are calculated considering the contri- 
bution from all the shocked region, and taking into account the 
relevant relativistic transformations of the radiation field, and 
the different photon arrival times to the different observers. 
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Fig. 2. — Light curves of model 3, for Oq = 0.2, E52 = no = 2 = 1, 
p = 2.5, ee = 0.1, es = 0.01, and u = 5 ■ lO" Hz. The inset shows 
the same light curves for model 2, where the same traces correspond to the 
same viewing angles 9chs- 
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Figure ^ shows the Ught curves of models 3, while the inset 
provides the light curves of model 2, for the same set of param- 
eters. In model 3, the peak of the light curves for 9ohs > Gq 
is flatter compared to model 2, and is obtained at a somewhat 
latter time. The rise before the peak is not as sharp as in mod- 
els 1 or 2, since in model 3 there is some material at the sides 
of the jet with a moderate Lorentz factor (Granot et al. 2001; 
Piran & Granot 2001). The emission from this slower material 
tends to dominate the observed flux at early times for observers 
at 9ohs > &o, resulting in a gentler rise before the peak. The 
light curves for 6'obs > ^'o peak at a later time compared to 
model 2, and the flux during the decay stage grows faster with 
9ohs, since in model 3 the curvature of the shock front is larger 
and the emission occurs within a shell of finite width, resulting 
in a larger photon arrival time, and implying that smaller radii 
contribute to a given observer time. The light-curves for model 
2 & 3 are quantitatively similar for 6'obs < ^'o- 

The main advantage of this model is a reliable and rigorous 
treatment of the jet dynamics, which provides insight on the 
behavior of the jet and the corresponding light curves. Its main 
drawback is the long computational time it requires. 

3. LINEAR POLARIZATION 

While the afterglow emission from a spherical outflow is ex- 
pected to exhibit little or no linear polarization, as the polariza- 
tion from the different parts of the afterglow image cancel out, 
a jetted outflow breaks the circular symmetry of the afterglow 
image and may exhibit a polarization of up to < 20% (Ghis- 
ellini & Lazatti 1999; Sari 1999). One might therefore expect 
an even larger polarization for an observer at 6'obs > do- 

An isotropic magnetic field configuration in the local rest 
frame will produce no linear polarization. However, as the mag- 
netic field is most likely produced at the shock itself, one might 
expect the magnetic field perpendicular (B±) and parallel (i3||) 
to the shock direction, to have different magnitudes (Gruzinov 
1999; Sari 1999). We calculate the linearpolarization for model 
2 following Ghisellini & Lazatti (1999)g and using their nota- 
tions. They assume the magnetic field is strictly in the plane of 
the shock (B = B^). 
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Fig. 3. — The linear polarization for model 2 for several viewing angle 
and for the same parameters as in Figure hi 

Figure || shows the polarization as a function of time for dif- 
ferent flobs in terms of Pqq. For Pgo < the polarization is 

^Our calculation is different from Ghisellini & Lazatti in that we include lateral 
results for O^hs < "o ar£ different. 



along the plane containing the line of sight and the jet axis, 
wile for Pgo > it is rotated by 90° (for (S^) < 2{B\\) this 
is reversed, e.g. Sari 1999). A more isotropic magnetic field 
configuration would result in a smaller degree of polarization, 
so the value of the polarization in Figure [l| 40%) may be 
viewed as a rough upper limit. For 0.3 <^ 6'obs/6'o <^ 1-1 the 
polarization vanishes and reappears rotated by 90° around tjet. 
This behavior may occur again at a later time, but the subse- 
quent polarization is very low. For 1.1 < 0obs/^o ^1-6 the 
polarization has two peaks, the first higher than the second. For 
(^ohs/Oo ^ 1-1 the polarization is largest near the peak in the 
light curve, and decreases quite slowly with time, while the 
peak polarization shows a very weak dependence on 6'obs, and 
is about a factor of 2 larger than for 6'obs — Oq. 

4. ORPHAN AFTERGLOWS 

If GRB jets have well defined edges, the prompt gamma- 
ray flux drops very sharply outside the opening of the jet, and 
the prompt burst will be very hard to detect from 6'obs > ^'o- 
On the other hand, the afterglow emission may be detected 
out to 9dct ~ a few 9q, where the exact value of f?dct de- 
pends on the jet parameters (including its redshift), the ob- 
served band and the limiting flux for detection. Jetted GRBs 
with 6*0 < 6'obs < ^'dct are expected to be orphan afterglows 
(i.e. detectable in the optical but not in gamma-rays). 

It has been argued by Dalai et al. (2002) that ^dot/^'o ~ 
const for 9o <C 1, so that the detection rate of orphan after- 
glows Norph (associated with off-axis jets) will be a constant 

[namely (6'dot/^o)^] times the GRB detection rate Nqrb' ™d 
thereby a comparison between these two rates will not constrain 
^0 or the true rate of GRBs Nqrb. This result was obtained as- 
suming a constant flux, F^{t jet), at tjct for 0obs = 0. However, 
afterglow observations suggest that the total energy in the jet, 
Ejct, is roughly constant (Frail et al. 2001, Panaitescu & Ku- 
mar 2001, Piran et al. 2001) while Fi,{tjct) varies over a wider 
range. In fact, for Ejct — const, simple jet models (R-SPH99) 
predict that the hydrodynamical evolution of the jet (and there- 
fore the light curves for all ^obs) becomes independent of 9q 
once the jet enters the phase of exponential lateral expansion 
with radius. This corresponds to t > tjct for 6'obs < Sq, and 
to t > (6'obs /^?o)^ijct for 9ohs > 00, which includes the time 
around or after the peak in the light curve. This implies that 
for Ejct — const, we have 6'dot = const = 0dct,o (rather than 
Odct/Oo ^ const) for da < 6'dct- For 6'o ^ 6'dct,o naturally 6'dot 
is larger 9o if the afterglow is detectable from 6'obs < ^'o, and the 
solid angle between da and 9dct, i^orph = cos 9a — cos 9dct ~ 
(^dot ^ dl)/2(m which we have detectable orphan afterglows) 
remains approximately constant. Furthermore, the distribution 
of 9q inferred from observations (Frail et al. 2001, Panaitescu 
& Kumar 2001) is sharply peaked at low 9a (~ 2 — 3°). This 
suggests that most of the orphan afterglows that would be de- 
tected should have 9q ^ 3°. For a reasonable limiting mag- 
nitude for detection, this implies (?dct ^ ^^o in most cases, 
and therefore ^dot ~ ^'det.o- For example, for model 2 with 
-Ejct = 5 ■ 10^° ergs (assuming a double sided jet), = 0.1, 
eB = 0.01, p = 2.5, na = z = \ and a limiting magnitude 
for detection of i? = 24 we obtain 12.4° < 6'dct < 23° and 
0.023 < l^orph/47r < 0.045 for 2° < 9o < 15°. If indeed 
i^orph ~ const, then Norph should provide a good estimate 

expansion of the jet and integration over equal photon arrival time surfaces, so our 
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of the true GRB rate, Ng^°B = (47r/norph)Norph- The aver- 
age beaming fraction = (1 — cos 6*0) « {Sq)/2 is given by 

h = NtAe/fe = (^^orph/47r)Nde^B/Ndet^. 

5. GRB 980425 /SN1998BW 

On April 25, 1998, a Gamma-Ray Burst was detected by 
Beppo SAX and CGRO. The burst consisted of a single wide 
peak of duration 30 s, peak flux in 24-1820 keV band of 3- 10"'^ 
erg cm~^ s"'^, and fluence of 4.4 • 10~^ erg cm~^ (Soffitta et al. 
1998, Kippen et al. 1998). The burst had no detectable emis- 
sion above 300 keV. The burst spectrum was a broken power- 
law with break at 148±33 keV, and the high energy power-law 
photon index of -3.8 ± 0.7 (see Galama et al. 1998). These 
values are not unusual for GRBs. 

A bright Type Ic supernova, SN 1998bw, located at z = 
0.0085, was detected within 8 arc minutes of GRB 980425. 
From the extrapolation of optical Ught curves Galama et al. 
(1998) suggested that the SN went off within a day of the 
GRB, thereby implying a possible connection between the two 
events. The probabihty of this association is strengthened by 
the uniquely peculiar Ught curve and spectrum of the SN (e.g., 
Patat et al 2001). Early on, Woosley et al. (1999) gave argu- 
ments why SN 1998bw might be a SN exploded by a jet and 
therefore possibly associated with a GRB. This would arise, for 
instance, in the coUapsar model (Woosley 1993; MacFadyen 
and Woosley 1999; Paczynski 1998). 

If indeed the two events are associated, then the total 
isotropic equivalent of energy release in 7-rays for GRB 
980425 is E^^^^o = 8.5 • 10^^ erg, or a factor of - 10^ smaller 
than the energy for an average cosmological GRB. This could 
explained in two ways. 

5.1. Sharp Edged, Homogeneous Jet Seen at dobs > ^0 

If GRB 980425 was viewed at 9ohs > &q it might explain its 
small E^^iso- For a GRB arising from a jet with 7 independent 
of 9 and sharp edges, the observed energy falls off rapidly for 
^^obs > ^0, in fact as &® where h = 7(0obs — ^^o)- Moreover, 
for an observer at Oohs > So the peak of the spectrum is lower 
by a factor 6^, and the burst duration longer by the same factor, 
compared to an observer at ^obs < ^0 • Therefore the low E^^iso 
of 980425 implies }fi - 10^ and Oohs ^ Oq + 5°(7/50)-i.' If 
GRB 980425 were at a cosmological distance and seen from 
^obs < Oq, the peak of the spectrum and the burst duration 
would be ~ 1 [3/ ( H- z)] MeV and 4 [( 1 -|- z) /3] s, respectively. 

A second constraint is set by the condition that the optical 
afterglow is dimmer than SN 1998bw, which had a luminos- 
ity of 2 X 10^^ erg at 1 day, rose to lO"'"' erg at 14 — 19 days 
(Galama et al. 1998), then decayed at 0.017 mag/day (Patat 
et al. 2001). To compare it with the afterglow luminosity, we 
shall use the typical properties of the afterglows whose optical 
light-curves exhibited breaks: i) average jet break tj^t ^ 0.5 
days in their rest- frame, ii) average optical luminosity flux 
LiOobs = 0,tjet) 2 X 10*^ erg at tjet, in) L oc t~" with 
a ^ 2 att > tjet, iv) L^, oc with /? 1, at optical fre- 
quencies. 

Using Model 1 described in section §2.1, it can be shown that 
for an observer at ^obs ^ 2^o the time and optical luminosity at 
the Ught-curve peak are 
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From these equations it can be shown that for ^obs ^ 3 the 
peak afterglow luminosity Lpeak{tpeak) is a factor 3 lower 
than L{t) of SN 1998bw. During the decay phase, the afterglow 
luminosity remains below that of the SN until about 400 days, 
when they become comparable. We note that Patat et al. (2001) 
report a flattening of the SN 1998bw decay after 375 day, to a 
dimming rate of ^ 0.009 mag/day, which they interpret as the 
settling in of the ^^Co decay or the existence of other sources. 

Together with the above constraint, ^obs ^ ^o + 5°(7/50)~^, 
imposed by the fluence of GRB 980425, the condition that the 
afterglow emission does not exceed that of the SN 1998bw 
leads to ^0 ^ 3°(7/50)-i and 3^o ^ ^obs ^ S°i^/50)-\ 

5.2. Inhomogeneous Jet Seen Off-Axis 

Another possibility is that the jet does not have sharp edges, 
but wings of lower energy and Lorentz factor that extend to 
large 9. Such a picture of the jet was suggested by Woosley 
et al. (1999) and is consistent with the relativistic studies of 
the collapsar model by Zhang, Woosley, & MacFadyen (2002). 
GRB 980425 would be then be produced by material with 
7 '--^ 10 moving in our direction. 

When seen on-axis, an afterglow with the above properties 
has L ~ 5 X 10^^ erg at 1 day and L - 2 x lO'*^ g^g at 16 
days, i.e. a "typical" afterglow would be 200 times brighter 
than SN 1998bw when its first observation was made and sev- 
eral times dimmer when SN 1998bw peaked. All other param- 
eters remaining the same, an afterglow emission dimmer than 
that of the SN 1998bw at 1 day requires an energy per solid 
angle ejot(^^obs) in the direction toward the observer satisfying 

^iOohsWlks - (200)-3/(P+3)^jet - 2 • 10^9 erg. As in the pre- 
vious subsection, the observer location satisfies ^obs ^ 39o so 
that the optical emission of the main jet of opening Oq is below 
the supernova light-curve. 

6. DISCUSSION 

We have presented the calculation of Ught-curves from a rel- 
ativistic jet for an arbitrary location of the observer; much of the 
work in this letter is for an observer located outside the initial 
jet opening, ^obs > ^o- We have considered three different jet 
models of increasing sophistication; the simplest being a point 
source moving along the jet axis (§2.1), and the most sophisti- 
cated is 2D hydrodynamical simulation (§2.3). The basic quali- 
tative features of the hght-curves are similar in all three models, 
for ^obs > ^0- Moreover, the uniform jet model (model 2, §2.2) 
is in rough quantitative agreement with the hydro-model. 

We find that "orphan" optical afterglows associated with off- 
axis jets can be observed up to a constant 9ohs, rather than a 
constant ^obs/^o as suggested by Dalai et al. (2002), if one 
assumes a constant energy in the jet, rather than a constant 
flux at the time of the jet break for an on-axis observer. This 
implies that future surveys for orphan afterglows may provide 
valuable data for the the distribution of jet opening angles 
and the true event rate of GRBs. The orphan optical events dis- 
cussed here can be identified from the initial rise during which 
the spectral slope is typically /3 > 0, followed by a decay, on 
a time scale of ~ 1 — 30 days, and may show a large degree 
of linear polarization (^ 40%). The detection of such orphan 
afterglows may provide a new line of evidence in favor of jet- 
ted outflows in GRBs. Recently Huang, Dai and Lu (2001) 
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have considered another scenario (failed GRBs) for producing 
orphan afterglows; this would increase the detection rate of or- 
phan afterglows. A good monitoring of optical transients may 
help distinguish failed GRBs from jets seen at 6'obs > 6*0^ 
improve our understanding of them. 
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